Phase Fluctuations in High Temperature Superconductors 
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Within the phase fluctuation picture for the pseudogap state of a high-Tc superconductor, we 
incorporate the phase fluctuations generated by the classical XY model with the Bogoliubov-de 
Gennes formalism utilizing a field-theoretical method. This picture delineates the inhomogeneous 
characteristics of local order parameters observed in high-Tc superconductors above Tc. We also 
compute the local density of states near a non-magnetic impurity with a strong scattering potential. 
The resonance peak smoothly evolves as temperature increases through Tc without showing any 
sudden broadening, which is consistent with recent experimental findings. 
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One of the most defining features of high-Tc super- 
conductors, the pseudogapstate, has drawn very intense 
attention in recent years p|. Various theoretical mod- 
els have attempted to explain this state; they gener- 
ally fall into two categories. One is thepre-formed pair 
model, or phase fluctuation scenario 
this picture, the Cooper pairs continue to exist above 
Tc, in the pseudogap state below T* without phase co- 
herence, where the superfluid stiffness is zero. A sec- 
ond category is the competing gap scenario, where, for 
example, a density-wave gap coexists with the super- 
conducting gap d, H [13, These distinctions can 
be fairly subtle, and the nature of this phase contin- 
ues to attract growing interest. Angle-resolved photo- 
emission spectroscopy experiments [1 2| have supported 
the phase fluctuation scenario since the pseudogap seems 
to evolve into the superconducting gap as the tempera- 
ture changes through T^. Further support can be found 
in the Nernst experiments [l3. 14 1, where a large Nernst 
signal induced by vortex motions is observed above T^ in 
hole-doped cuprates. This result has been interpreted 
as evidence for phase fluctuations. In the pseudogap 
state, from the viewpoint of the phase fluctuation sce- 
nario, one can conceive a spatially dependent order pa- 
rameter A(r) = |A|e'*'^''^ at a location r 15|. In fact, 
an intrinsic inhomogeneitv [l^ [l7j of local order param- 
eters has been observed even at a low temperature. In a 
two-dimensional theory 0, 0, H, d, 0|, the transition into 
the superconducting state at T^ is not BCS-like. Instead, 
it is the Kosterlitz-Thouless (KT) transition Con- 
sequently, we identify Tkt — Tc, which depends on the 
superfluid stiffness. Denoting Tmf as the temperature 
when Cooper pairs preform, we also define Tmf = T*, 
which is determined by the strength of the pairing po- 
tential. 

Recently, spatial variations of the ord er p arameter 
were visualized in topographic images [3, [loj. Gomes 
et al. [l9| have examined statistically the evolution of the 
order parameter in atomic scale of Bi2Sr2CaCu208+a: for 
various doping cases with increasing T from a low T 



through Tc to a high T. They illustrated the distribu- 
tion D as well as the probability P of local order param- 
eters. The distribution D{A) is an ordered list of local 
order parameters, which describes how many sites have 
the order parameter less than A. Thus D{A -t- e) — -D(A) 
gives the number of sites with the order parameter be- 
tween A -I- e and A. Since P(A) is the probability to 
find a local order parameter of A, P(A) oc di^D{uj)\A- 
On the other hand, Chatterjee et aZ.|20|] observed that 
the impurity resonance survives for T > Tc in the 
study of the scanning tunneling microscop y ( STM) for 
Bi2_yPbySr2CaCu208+^ . The STM studyjlj for high- 
Tc superconductors is performed usually at a low T <C Tc. 
Since Tc of this compound is 15K, they carried out the 
T dependence of the impurity resonance through Tc. In 
this paper, we wish to demonstrate that these recent 
findings 19|,l20| are understandable, at least qualitatively. 



within the framework of the phase fluctuation scenario. 

In the theoretical formulation, we separate the order 
parameter into the mean-field value and the phase fluc- 
tuation part using the fleld-theoretical approach. The 
mean-field value is determined by the the Bogoliubov-de 
Gennes(BdG) formalism while phase fluctuations will be 
generated by the classical XY model. Let us start with 
the Hamiltonian of a d-wave superconductor with a def- 



inition of ^pj = (Cj|, Cii) and the Pauli matrices fa and 

f± = (fi ± iT2)/2, 

- J2^^,J4'lr+^PJi'}f-i>, (1) 



where tij is the hopping amplitude, Ui is an impurity 
potential at a site i, is the chemical potential, Vij is 
a pairing potential, and a = —a with a =\, |. The 
symbol means a sum over nearest neighbor pairs. 

Since no magnetic field is included, t^j = i^i+^^i, where 
i = ±x, ±y are the unit vectors along the x{y) direc- 
tion. Using a set of impurity sites I, one can write Ui 
as UiSi^i For the nearest neighbor pairing, Vij — 
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VSi+s,j- We set the lattice constant to be unity and use 
units such that h = ks = 1- The partition function of Ti. 
in the path integral is Z = J V^Vip'^ ^-■31^,4'^] where 

= IdT[J2^^^l^r^'^+H\ and Vij = U^dA■ 
The range for the integral over imaginary time r in the 
action S is from to f3 ~ 1/T. 



Introducing the Hubbard-Stratonovich transformation 
with auxiliary field (jjij as a field-theoretical method to 
deal with phase fluctuations, we obtain the partition 
function as follows: Z = / I?0P0tp^x>^t g-S[0,</.t,V,Vt]^ 

where, for the static case: dripi — 0, 



dr 



<«j> 



(2) 



r 



In the saddle point approximation, SS[(j), (f)^ , ip, ^p^/S(j)^ = 
0, the auxiliary field can be identified as the order 
parameter; namely, (pij — Vijtp^jT-'ipi = Ay. Now 

the BdG Hamiltonian is 

<i,j> i 

- E {^^oi^l^+^o + h.c) . (3) 
The effective partition function can be written as 



where 7Y, 
ing order 

(W4)E„ 



parameter is evaluated 



n n* I „,n*..n 



tanh(i;„/2T). 



The bond- 
as A,;^ 

lo- 



The 



cal d-wave order parameter at a site i is given by 



A'* = i y - A 

ously speaking. Ay 
Consequently, A^ = 



A,, 



where Aj^ = |A.y|e^^'^. 



exp 



ZsdG (4) 



where Z 



BdG 



= Tr[e 



-I3He 



'] is the partition function 



corresponding to the BdG Hamiltonian. Note that ZBdo 



\Aij I exp \i{9ij + 9ij)\ . The phase 
Oij associated with the mean-field value is determined 
within the BdG calculations. On the other hand, Oij 
is generated based on the XY model as mentioned. 
Assuming the effective phase of the local d-wave order 
parameter varies slowly over a distance of the coherence 
length^, it is reasonable to define Oij = [Oi + Oj] /2, 
which corresponds to the phase of the center of mass of 
a Cooper pair in the continuum limit. 

Let us now calculate the average of (Ay) over phase 
fluctuations. It is straightforward to show that (Ay) = 



depends on A^^ — |Ay |e*^'^, where Oij represents phase ^iji^'' '^): where 



fluctuations (see below). However, as a mean field ap- 
proximation, the BdG formalism does not consider phase 
fluctuations. To incorporate the fluctuations into the 
framework, we utilize the two-dimensional XY model|22 
with HxY = ->/xy I]<iy> cos(6'i 

the coupling strength (or the superfluid stiffness) and Oi 
is the angle made by a classical spin at a site i. Following 
Ref . [2^1 ■ we perform simulations using the Monte Carlo 
method with the Metropolis algorithm. The lattice size 
in our calculations is 24 x 24 with periodic boundary con- 
ditions. The minimum length of phase fluctuations is set 
to be 2. Starting from a random configuration, the XY 
model system is cooled down to a working T. It is known 
that Ttk ~ 0.9 JxY in numerical simulations. To bench- 
mark our simulations, we also compute the specific heat 
of the system, which shows similar behavior to the one 
computed in Ref. [2^ . 

The mean-field value Aij of the order parameter is 
obtained self-consistently by solving the BdG equation: 



/ VO e'^^^ZBdcjA, 

J VOZBdG[^,0] 



(6) 
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A* 



A,. 

-n* 



En 



,), where Jxy is with ZedG = n„(l + e"^^")- As in Refs.[i, 0], 
(e*^*^) is numerically evaluated using the Monte Carlo 
method. Generating a phase configuration {0^} at 
T and diagonalizing the BdG Hamiltonian with Ay , 
we obtain £"„ corresponding to the given configura- 
tion. In the averaging process, the total number of 
phase configurations of an ensemble at T is about 10''. 
The local d-wave order parameter now is evaluated as 
(Aj ) = i '}Zj{Aij) \6j^i+x + 5j^i-x - 5j,i+y - 5j^i-y]. Fol- 
lowing this procedure yields a statistical distribution D of 
local order parameters A = |(Af)|. As explained earlier, 
the probability P{A) is given by d^D{ijj)\i:^. In Fig. 1 
we plot P{A) (in an arbitrary unit) at various T in the 
units of the hopping amplitude. Other parameters used 
are V = 2.05, /i = 0, and Jxy = 0.11 and impurities 
are not considered. At low T, P(A) shows a sharp peak 
at A = 0.25, which is the mean-field value at T = 0, 
because phase fluctuations are not significant. Even near 
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FIG. 1: (Color online) Probability P(A) (in an arbitrary unit) 
to find the value of local order parameters A at various T in 
the unit of the hopping amplitude. As T increases, the prob- 
ability distribution spreads. However, P(A) shows a sharp 
distribution again when T Tmf- The inset shows that 
the spatial average (symbol) of the local order parameters 
diminishes as T increases compared with the mean-field gap 
(curve). 



Tc ~ jTmf, -P(A) is confined within A ~ 0.21 - 0.24. 
As T increases above Tc, the probability spreads; how- 
ever, it sharpens back when T — > T]\ip because the mean- 
field gap diminishes rapidly. Since the median of P(A) 
moves towards zero, the spatial average of the local order 
parameters decreases. The inset of Fig. 1 demonstrates 
such a behavior compared with its mean-field counter- 
part. In particular, the spatial average does not change 
much until T > 1.5Tc. This is consistent with a vari- 
ation of the density of states computed in Ref.Q]- We 
also performed the same calculations choosing V = 0.92 
for which Tc ~ ^Tmf and obtained qualitatively similar 
results. Gomes et al.^^ measured local c?-wave order pa- 
rameters of Bi2Sr2CaCu208-i-a; to show the order param- 
eter distribution. They also presented the probability for 
the overdoped Bi2Sr2CaCu20g+2; in a histogram with an 
interval of 2r7iey, which is equivalent to e in our defini- 
tions. The histogram illustrates the probability spreads 
as T increases while it becomes a sharp peak at a high 
T. It is also indicated that the median of the probabil- 
ity moves towards zero with increasing T. The detailed 
behavior of the probability in Ref. [l^l is not completely 
identical to that shown in Fig. 1. Nonetheless, impor- 
tant characteristics of the histogram are congruent with 
Fig. 1. 

Another interesting experimental observation usin g th e 
STM has been recently reported by Chatterjee et aZ.[20l|. 
They observed the evolution of the impurity resonance 
as T increases from low T to high T(> Tc) and found 
that the resonance survives above Tc- Moreover, the res- 



onance evolves smoothly with increasing T and does not 
show any sudden broadening near Tc- In fact, it was ar- 
gued that such a sudden broadening near Tc would occur 
in the phase fluctuation scenario j23|. This led Chatter- 
jee et al. to conclude that the phase fluctuation sce- 
nario is not consistent with their findings. However, in 
Ref. [2^ the phase of the local order parameters is as- 
sumed to vary on the length scale of the London pen- 
etration depth so that the quasiclassical approximation 
is applicable. Within our framework, such an assump- 
tion corresponds to weak phase fluctuations and, thus, Tc 
would not be much less than T* . The only assumption we 
made following Refs. [1, 01 is that the phase varies slowly 
over the coherence length without any larger length scale 
employed, and the BdG formalism is applied. Conse- 
quently, the phase varies at any length scale greater than 
the lattice constant, and any degree of fluctuations can 
be considered. 

The local differential conductance is described by the 
local density of states (LDOS) in the theoretical calcula- 
tions. We apply the supercell technique to compute the 
LDOS as in Refs. 24, 25]. If the lattice size is x Ny 
as a unit cell and the number of unit cells is x My, 
then the total size becomes N^M^ x NyMy. We also 
introduce the quasimomenta Px{y) = 1v~^^^mY7' where 
^x(y) = 0, 1, • • • Mx(^y) — 1. It is obvious that Zbclg, now, 
depends on Px{y)- In the presence of phase fluctuations, 
one can think of two different supercell techniques de- 
pending on how to associate the fluctuations with Px(y)- 
One way is to start with the same random conflguration 
of phases regardless of Px{y) while the other would be 
to use different random configurations for different val- 
ues of Px(y)- Nonetheless between the two methods we 
found no practical difference in the LDOS averaged over 
the ensemble. In particular, for T > Tc we obtain iden- 
tical results. We obtain the averaged LDOS as follows: 
the LDOS is calculated for a given phase configuration 
and its average is computed over the ensemble. For a 
single impurity, we choose a strong impurity potential 
U = 100, which is close to the unitary limit [2q. Other 
parameters are V — 1.1, ^ — 0, and Jxy = 0.033. For 
the lattice size, Nx = Ny = 24 and Mx = My = 10. In 
Fig. 2, we plot the LDOS at the nearest neighbor site to 
the impurity (solid curves) as well as at a site far away 
from the impurity (dashed curves) for T = 0.01, 0.02, 
0.03 and 0.04. Note that the impurity resonances sur- 
vive above Tc ~ 0.03. Moreover, the resonance does not 
change much with increasing T near Tc- The displayed 
behavior is consistent with the experimental findings [20|. 
It is illustrated further in Fig. 3 for = 1.1 and for a 
weaker pairing potential V — 0.6 at T = 0.02, 0.03, and 
0.04. This indicates that the impurity resonance above 
Tc does not strongly depend on parameters. From the 
theoretical point of view, it is originated from the fact 
that for these cases phase fluctuations (6*^) cannot dis- 
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FIG. 2: (Color online) LDOS at the nearest neighbor of the 
impurity (solid curves) and far away from it (dashed curves) 
for various T. Regardless of the value of V , the impurity 
resonance evolves smoothly with increasing T and survives 
above Tc ~ 0.03. The relevant parameters are V = 1.1, U = 
100, M = 0, and Jxy = 0.033. 



FIG. 3: (Color online) LDOS at the nearest neighbor site of 
the impurity near Tc for V = 1.1 and 0.6. Other parameters 
are same as in Fig. 2. The resonance smoothly evolves with 
increasing T from 0.02 to 0.04 through Tc. 



turb significantly the phase difference of the mean-field 
value 7r(= 9i^i+y - 9i,i+x)- 

In conclusion, considering recent experiments on the 
imaging of the inhomogeneous local order parameters 
above Tc as well as the impurity resonance peak in high- 
Tc superconductors, we calculated a distribution of the 
local order parameters and the local density of states near 
an impurity based on the phase fiuctuation scenario. The 
order parameter distribution is qualitatively agreed with 
experimental observations. Randomly distributed impu- 
rities may be necessary for more realistic model calcu- 
lations. A smooth evolution of the impurity resonance 
peak, through Tc, seen in the STM measurement also 
supports this picture. The resonance peak computed 
within our framework does not show any sudden broaden- 
ing near Tc and consistent with the experimental results. 

We acknowledge N.P. Ong for useful discussions. This 
work was supported by the Robert A. Welch Foundation 
and the Texas Center for Superconductivity at the Uni- 
versity of Houston through the State of Texas, the Min- 
istry of Science and Technology of China (973 project No: 
2009CB929200), NSFC Grant No. 10874032, and the 
Program for Professor of Special Appointment (Eastern 
Scholar) at Shanghai Institutions of Higher Learning. 



[1] T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 
[2] V.J. Emery and S.A. Kivelson, Nature 374, 434 (1995). 
[3] S.G. Sharapov, H. Beck, and V.M. Loktev, Phys. Rev. B 

64, 134519 (2001). 
[4] T. Eckl, D.J. Scalapino, E. Arrigoni, and W. Hanke, 



Phys. Rev. B 66, 140510 (2002). 
[5] W. Kim and J. Carbotte, Europhys. lett. 59, 761 (2002). 
[6] L. Benfatto, A. Toschi, and S. Caprara, Phys. Rev. B 69, 

184510 (2004). 

[7] Y. Dubi, Y. Meir, and Y. Avishai, Nature 449 876 (2007). 
[8] S. Chakravarty, R.B. Laughlin, D.K. Morr, and C. 

Nayak, Phys. Rev. B 63, 094503 (2001). 
[9] J.X. Zhu, W. Kim, C.S. Ting, and J.P. Carbotte, Phys. 
Rev. Lett. 87 197001 (2001). 
[10] S. Sachdev, Rev. Mod. Phys. 75 913 (2003). 
[11] S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75 1201 (2003). 
[12] H. Ding, J. R. Engelbrecht, Z. Wang, J. C. Campuzano, 
S.-C. Wang, H.-B. Yang, R. Rogan, T. Takahashi, K. 
Kadowaki, and D. G. Hinks, Phys. Rev. Lett. 87, 227001 

(2001) . 

[13] Z.A. Xu, N.P. Ong, Y. Wang, T. Kakeshita, S. Uchida, 
Nature 406, 486 (2000). 

[14] Y. Wang, N.P. Ong, Z.A. Xu, T. Kakeshita, S. Uchida, 
D.A. Bonn, R. Liang, and W.N. Hardy, Phys. Rev. Lett. 
95, 247002 (2005). 

[15] The thermal average of the order parameter depends on 
phase fluctuations at a given temperature (T): (A(r)) = 
|A|(e'^'""'). Since (e'^'*''^) ~ e'^^^''^ for a weak fluctua- 
tion, |{A(r))| ~ |A|. On the other hand, |(A(r))| vanishes 
when the fluctuation is strong because (e**'"'') 0. 

[16] S.H. Pan, J.P. O'Neal, R.L. Badzey, C. Chamon, H. Ding, 
J.R. Engelbrecht, Z. Wang, H. Eisaki, S. Uchida, A.K. 
Gupta, K.-W. Ng, E.W. Hudson, K.M. Lang, and J.C. 
Davis, Nature 413, 282 (2001). 

[17] K.M. Lang, V. Madhavan, J.E. Hoffman, E.W. Hudson, 
H. Eisaki, S. Uchida, and J.C. Davis, Nature 415, 412 

(2002) . 

[18] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 
(1973). 

[19] K.K. Gomes, A.N. Pasupathy, A. Pushp, S. Ono, Y. 
Ando, and A. Yazdani, Nature 447, 569 (2007). 



5 



[20] K. Chatterjeo, M.C. Boyor, W.D. Wise, T. Kondo, T. 
Takcuchi, H. Ikuta, and E.W. Hudson, Nature Physics 
4, 108 (2008). 

[21] S.H. Pan, E.W. Hudson, K.M. Lang, H. Eisaki, S. 

Uchida, and J.C. Davis, Nature 403, 746 (2000). 
[22] J. Tobochnik and G.V. Chester, Phys. Rev. B 20, 3761 

(1979). 

[23] Q.H. Wang, Phys. Rev. Lett. 88 057002 (2002). 



[24] J.X. Zhu, C.S. Ting, and C.R. Hu, Phys. Rev. B 62, 6027 

(2000). 

[25] Y. Chen and C.S. Ting, Phys. Rev. Lett. 92, 077203 
(2004). 

[26] M.L Salkola, A. V. Balatsky, and D.J. Scalapino, Phys. 
Rev. Lett. 77, 1841 (1996). 



